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Ph; abstract 

o , 

^ . Image modeling and simulation are critical to extending the limits of leading edge lithography technologies used 

O ' for IC making. Simultaneous source mask optimization (SMO) has become an important objective in the field of 

' ^ . computational lithography. SMO is considered essential to extending immersion lithography beyond the 45nm 
node. However, SMO is computationally extremely challenging and time-consuming. The key challenges are due 

*-£h , to run time vs. accuracy tradeoffs of the imaging models used for the computational lithography. 

'— 1 , We present a new technique to be incorporated in the SMO flow. This new approach is based on the reduced 

basis method (RBM) applied to the simulation of light transmission through the lithography masks. It provides a 
, rigorous approximation to the exact lithographical problem, based on fully vectorial Maxwell's equations. Using 
(S) ' the reduced basis method, the optimization process is divided into an offline and an online steps. In the offline 
, step, a RBM model with variable geometrical parameters is built self-adaptivcly and using a Finite Element 
(FEM) based solver. In the online step, the RBM model can be solved very fast for arbitrary illumination 
and geometrical parameters, such as dimensions of OPC features, line widths, etc. This approach dramatically 
reduces computational costs of the optimization procedure while providing accuracy superior to the approaches 
' ^ , involving simplified mask models. RBM furthermore provides rigorous error estimators, which assure the quality 
' and reliability of the reduced basis solutions. 

We apply the reduced basis method to a 3D SMO example. We quantify performance, computational costs 
• i-H . and accuracy of our method. 
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1. INTRODUCTION 



The importance of numerical simulations in the production process of integrated curcuits has grown rapidly 
over the past years. Accuracy requirements for lithography and scatteromctry simulations thereby increase 
hand in hand with decreasing feature sizes and fabrication tolerances. In many application areas like design and 
optimization of photomasks or scatteromctry, rigorous electromagnetic near- field computations have to be carried 
out for quantitatively correct results. With mask structures much smaller than the wavelength of incident light, 
approximative models like Kirchhoff do not give the correct physical solution. Usage of rigorous electromagnetic 
field (EMF) solvers, e.g., based on Finite Elements (FEM), Finite Difference Time Domain Methods (FDTD), 
or Rigorous Coupled Wave Analysis (RCWA), however comes at the price of higher computational costs. 

When solving inverse problems for scatterometric parameter reconstruction or numerical determination of an 
optimal mask layout, not only a single simulation has to be carried out, but multiple EMF solutions have to be 
computed. Then high computational times for a single forward solutions can make the complete simulation task 
infcasible. 

However, above applications have a specific feature. Usually the same basic layout has to be simulated 
multiple times for different values of geometrical parameters, e.g., line width or height, absorber edge angle, etc. 
The reduced basis method 1-5 can be applied to this setup. The solution process of the parametrized system 
is decomposed into an offline step which has to be performed only once and a fast online step. In the offline 
step the reduced basis is built self-adaptively by solving the underlying model rigorously several times. The full 
model is then projected onto the reduced basis. In the online step the assembled reduced system can be solved 
in the order of seconds independent on the size of the original problem. Furthermore, methods from the well 
established field of a posteriori error estimation of Finite Element theory 6 can be applied to the reduced basis 
method to assure reliability of the computed output and also for construction of the reduced basis. 7,8 

In the following we start with the description of an exemplary source mask optimization problem. We then 
give the key ideas and features of the reduced basis method, which is used for efficient solution of the SMO project 
in the final section. There we also perform a convergence analysis of our method quantifying its performance. 

2. SOURCE MASK OPTIMIZATION 

In this section we describe a source mask optimization problem, which we want to solve with the reduced basis 
method. As explained in the abstract, source mask optimization is an important objective for realization of ICs 
of present and future technology nodes. In SMO the shape of the photomask is not optimized for a fixed source 
like, e.g., in conventional optical proximity correction, 9 but the source itself is variable and included into the 
optimization process. The goal is to find an optimal source together with the optimal mask simultaneously. The 
combined degrees of freedom of source and mask enlarge the dimension of the parameter space and therewith 
also the computational costs for determination of the optimal parameter values significantly. Furthermore, since 
structures on the mask are getting much smaller than the wavelength of incident light, approximative models 
like Kirchhoff are not valid and rigorous solutions have to be computed during optimization. SMO is therefore 
considered computationally extremely challenging and time consuming. 

The optimization problem we want to tackle is an exemplary cutting lithography step, 10 used for hole 
patterning of a line-and-space pattern, see Fig. 1(a). Our goal is maximization of the process window. Parameters 
of the system are given in Table 1. 

SMO for this setup, based on Kirchhoff 's method has been presented in. 10 

The photomask and parametric illumination of this example are shown in Fig. 1(b), (c). For determination of 
the process window, given exemplary in Fig. 1(d), we measure width and height of the holes created on the wafer. 
Now we parametrize the mask and illumination. The mask topology is described by three geometrical parameters 
visualized in Fig. 2. As variable source we introduce four parameters to the quadrupole illumination, depicted 
in Fig. 3. Frccform optimization of the source leads to no principal changes of the reduced basis flow, however 
in this demonstrative study we only consider a parametric source. Each point in the pixclatcd representation of 
the illumination stands for two plane waves with orthogonal polarization to model unpolarized light. 



(a) 



(b) 




Figure 1. 

(a) Target cutting hole structure on wafer; dashed line denotes unit cell of layout; arrows denote position for 
determination of hole width and height for process window definition, (b) Top image of a unit cell of the periodic 
photomask, (c) Pixelated quadrupole illumination, (d) Sample combined process window for target width and 
height of cutting hole (±10%). Vertical black line: height of the process window at zero defocus; horizontal black 
line: width of the process window at midpoint of vertical line. For optimization, the size of the process window 
was taken as the product of the length of the sketched black lines. 



parameter value 
wavelength A 
NA 

magnification 

""substrate 
^substrate 
"absorber 
^absorber 

Table 1. Parameters of source mask optimization example, 
(a) (b) (c) (d) 




Figure 2. Geometrically parametrized mask, (a) Basic layout, (b) Variable hole width w. (c) Variable hole height h. 
(d) Variable assist feature d. (c.f. Fig. 1(b)) 

Our optimization task is the maximization of the process window in the seven-dimensional parameter space. 
For determination of a single process window, we first compute the electromagnetic near-field for fixed source and 
mask geometry. From this we deduce far field coefficients, aerial image, and contour plots on the wafer. From 
these we compute the process window. Here, computation of the near-field is by far the most expensive step. 
Discretization with Finite Elements gives a system with Af = 341440 unknowns, which results in single CPU 
computational times of the order of 10, 000 sec for 18 plane wave sources, consisting of nine incoming directions 
times two orthogonal polarizations, see Fig. 3(b). All these incoming directions have same Bloch-periodicity, i.e., 
phase jump from left to right and top to bottom boundary of the computational domain, and can be computed 
simultaneously in a Finite Element simulation. These sources are used in a local/hybrid 11 Hopkins approach to 
obtain the far field coefficients for the continuous (numerically pixelated) off-axis illumination. Thereby, the far 
field of an arbitrary incoming field is determined from the Finite Element solution which has the closest incoming 
direction. 

During optimization, several hundred of these forward solutions have to be computed, which makes the task 
basically infeasible, due to long run-times. 

In the following we describe the reduced basis method which will dramatically decrease the computational 
time for the optimization process. It replaces the near-field solution by a quasi-rigorous solution of controlled 
quality, which is obtained from a reduced model. We will give the key features and ideas of the reduced basis 
method in the following section. A detailed mathematical discussion can be found in. 

3. REDUCED BASIS METHOD 
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For rigorous numerical solution of Maxwell's equations, a discretization scheme like FEM, FDTD, or RCWA has 
to be applied. This results in a system of algebraic equations for the degrees of freedom which represent the 
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Figure 3. (a) Quadrupole illumination with four variable parameters; 61, 62: segment opening angles, ri, T2 inner segment 
radii, (b) Incoming directions in the cr-plane. Crosses denote incoming directions with equal Bloch-periodicity which are 
computed simultaneously in the Finite Element and reduced basis simulation. 

electromagnetic field. Especially in 3D, the computational times for determining the solution to the discretized 
equations can get very large. For inverse or optimization problems a large number of such solutions have to be 
computed for different parameter values of a parametrized system. This makes the solution of such real-time or 
many-query problems extremely time-consuming, inconvinicnt for a user, and often infeasible. 

The idea of the reduced basis method is the usage of precomputed solutions in an actual application. In 
a relatively expensive offline step a reduced model is constructed by multiple solution of the original problem 
for different parameter values. The assembled reduced model can then be solved very fast in the online phase. 
During an actual application, only the online step is evoked. At first glance, above strategy seems comparable to 
a table based ansatz where solutions arc computed for different parameter values offline followed by interpolation 
online. However, in the following we will see that the reduced basis method has great advantages regarding offline 
computational times and reliability of the solution. 

For application of the reduced basis method, a discretization with Finite Elements is the method of choice 
and builds its foundation. With Finite Elements, Maxwell's equations are transformed into a sparse matrix 
equation: 

Au = b, 

where A is the sparse system matrix, u the solution vector, and b the right hand side, which includes source fields, 
we call this system our "truth approximation" . For real- world problems the dimension AT of this system can 
be up to several millions, which leads to large computation times. Now consider a geometrically parametrized 
problem as depicted in Fig. 2. Since the domain of interest is parameter-dependent, the Finite Element system 
and the solution u will also become parameter- (y-) dependent: 

A v u v = b v . (1) 

During optimization we want to solve this equation for a large number of values for i>, which are elements of a 
parameter space D. 

The basic idea of the reduced basis method is depicted in Fig. 4. Suppose we have only one parameter which 
is varied. Then the solution it„ will stay on a (maybe complicated) one-dimensional sub-manifold hA 1 ^ in the 
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Figure 4. Sub-manifold of possible solutions M in Finite Element space X for a single parameter v. Three snapshot 
solutions u(yi) are depicted. 

high-dimensional Finite Element space. Suppose we would know a subspace U (AT x TV-matrix) of the Finite 
Element space, which gives a good approximation to the sub-manifold of all possible solutions: 

Then we can represent the solution u v of (1) by a linear combination of elements of U : 

u v « U\ u = u u , (2) 

where \ v is a N dimensional coefficient vector. Now we insert ansatz (2) into the discretized Maxwell's equations 
(1) and project the resulting system onto U. This gives: 

[U H A V U} K = u H K, (3) 

where the superscript H denotes the Hermitian transpose. This system of equations only has N unknowns, the 
coefficient A„, and can be solved very fast if N stays moderate 0(100). Equation (3) is the reduced basis system, 
U the reduced basis, and u u the reduced basis solution. An important issue is the construction of the reduced 
basis U, which we explain in the following section. Before that, we will comment on the features of the reduced 
basis method and compare its properties to table based approaches. 

In the online step of a table based ansatz, the physics of the underlying system is completely neglected. 
The output of interest for a specific parameter is simply computed by interpolation from output values of other 
parameters. In the reduced basis setting Eq. (3) can be interpreted as a scattering problem stated on the space 
of global ansatz functions U instead of the Finite Element space. 7 This means in the online step we still solve 
Maxwell's equations rigorously, but only on a reduced solution space. As we will see, this leads to a very fast 
exponential convergence of the method with increasing dimension of U and is a way to overcome the "curse of 
dimensionality" for high-dimensional parameter spaces which is a serious problem for all table based approaches. 

Another great advantage of the reduced basis method is the possibility to estimate the error of a reduced basis 
solution, without computing the exact Finite Element solution. This is the field of a posteriori error estimation, 
well investigated in the field of Finite Elements. The basic idea is an evaluation of the residuum of a reduced 
basis solution, which is defined by: 



From discretized Maxwell's equations (1) we see that the residuum vanishes for the exact solution u u . If wc 
insert a reduced basis solution into (4) and compute its residuum we can utilize ||fv|| as a measure for the error 
of the solution. Wc will demonstrate in our numerical example the good performance of this estimator. Wc 
can use this estimator to scan the parameter domain and identify regions, where the reduced system still needs 
improvement. Comparable methods do not exist for table based approaches. 

The error estimator is also used for self adaptive construction of the reduced basis approximation. In an 
application the user only specifies the parameter space D and the error estimator decides how to construct U as 
will be explained in the following. 

3.1. Construction of the reduced basis space 

The reduced basis space U has to be constructed such that it gives a good approximation to the space of all 
possible solutions M.^ of the parameter domain. In order to explain the construction algorithm, we first define 
a sequence of hierarchical subsets of the parameter domain D. Let i>i G D , i = 1, . . . , -ZV max . Then we define: 

Si ={ui,...,Vi} , i = l,...,N xaax . (5) 

These sets have the following property: 

Si C S 2 C ••• C S N C ••• C 5jv max . (6) 

The reduced basis space Un of dimension N is then defined by: 

Um =span{u(f) is a solution to (1) | v e Sn} , (7) 

hence it is spanned by solutions to the truth approximation for fixed parameter values. These solutions are called 
snapshot solutions, see Fig. 4. It is also possible to include first and higher derivatives of the field u with respect 
to parameters into the reduced basis space, which leads to so called Taylor and Hermite spaces. Here we notice 
the expensive "offline" costs of the reduced basis method. For construction of a reduced basis space of dimension 
N the truth approximation, i.e., the full problem, has to be solved N times. However, the reduced basis space 
only has to be assembled once and furthermore a parallclization of this process is possible. 

The final question is how to choose the snapshot parameters V{ in (5). The following greedy algorithm 3 results 
in reduced basis spaces of small dimension and very good approximation quality. 

First we define a training set Strain £ D of possible snapshot candidates. From this space we want to chose 
a number N of snapshot parameters for construction of the reduced basis. The first parameter v\ is chosen 
randomly. Then we construct a one dimensional reduced basis approximation corresponding to the snapshot 
u(y\). Now we evaluate the error estimator for this one dimensional reduced basis approximation on all candidate 
snapshots in the training sample 5 tra in and include the parameter value with the maximum error into the reduced 
basis because this is supposed to add a maximum of "new information" into the reduced basis. Then wc have 
a two dimensional reduced basis and the process is continued itcrativcly. The process can be stopped, e.g., if a 
certain maximum dimension is reached or if the error estimator gives sufficiently small bounds over the training 
set. 

4. NUMERICAL RESULTS 

In this section we perform the source mask optimization described in Sec. 2 with a reduced basis approximation 
for the electromagnetic near-field computation. First we specify the geometrical parameters of the mask, whose 
unit cell is shown in fig. 1(b). The pitch values are given by p x = p y = 704 nm. The geometrical parameter 
space D of the mask was chosen as: 



w E [120; 200] nm, 
h e [472; 632] nm, 
d E [24; 120] nm, 



(8) 
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Figure 5. Convergence of reduced basis solutions over a random parameter ensemble H; minimum, mean, and maximum 
errors are shown, (a) Error in i/(curl)-norm of near-field solution, (b) error of diffraction mode, (c) Scatter plot of 
estimated errors in dependence on true //(curl)-error over the random parameter ensemble and reduced basis dimensions 
f to 40. 



c.f. Fig. 2. 

For above parameter space, a reduced basis approximation of dimension 40 was constructed with the described 
greedy algorithm. Before we look at the optimization results we want to quantify the accuracy of the reduced 
basis approximation. Therefore, we randomly choose 50 parameters in the parameter domain (8) and compute the 
corresponding exact Finite Element solutions. This defines a random ensemble 5. Then we compare each of these 
exact solutions with the reduced basis results for the same parameters and increasing dimension of the reduced 
basis approximation: in Fig. 5(a) the errors of the reduced basis near-field approximations in H (curl , f2)-norm 
are shown, Fig. 5(b) gives the error of a reduced basis diffraction order, which is used for computation of the 
aerial image. The mean, minimum and maximum errors over the random ensemble S are shown. We observe 
exponential convergence of the errors with increasing reduced basis dimension. This means the more snapshots 
we include into the reduced basis system, the more accurate results we obtain. Reduced basis approximation 
converges to the exact Finite Element solution. The mean error for the output of interest, i.e., the diffraction 
orders is below 10 -4 for a dimension of only 40 which guarantees reliable results. 

As explained before, the high convergence rate is a major advantage of the reduced basis method. Only few 
snapshot solutions in the offline phase are necessary to obtain results with very small errors in the online phase. 
Next we look at the performance of the error estimator. Figure 5(c) shows a scatter plot of the estimated error in 
dependence on the true error over the random parameter ensemble. Thereby, the error of the near-field solutions 
given in 5(a) is given in dependence on the estimated error for each of the 50 random parameters and reduced 
basis dimensions 1 to 40. We observe excellent correlation of both quantities over several orders of magnitude 
of the true error. This assures reliability of the reduced basis approximation, since the error estimator can be 
utilized to check if the reduced basis solution has sufficient quality in the whole parameter space. Solution of the 
reduced basis system of dimension 40 takes approximately Is. Compared to a full Finite Element simulation, 
this results in a speed-up factor of 10,000. 

Due to the short online solution time of the reduced basis model, it is now very convenient to perform 
parameter scans and source mask optimizations. 

4.1. Source mask optimization 

As optimization routine we use a combination of global and local minimization algorithms. First we perform 
a scan with fixed source over the geometrical parameter domain (with 9-9-6 = 486 forward solutions) and 
determine mask geometries with large process windows. Figure 6 shows a poor (a,c,e) and the best (b,d,f) result 
for this run. We nicely observe a narrow waist of the hole pattern in the aerial image in Fig. 6(c) for the poor 
parameters, which is mainly due to the large value of d, c.f., Fig. 1(b). The target hole height and width can 
not be printed simultaneously for these parameters. 

For the global geometrical scan, we fixed the source to n = r 2 = 0.8 and b\ = b 2 = 45°. The "good" 
geometry candidates now serve as initial guesses for local optimization routines (we use a simplex algorithm, 
with a maximum number of iterations of 100), which are evoked on the common parameter domain of source 
and mask. 

The optimal process window is obtained for the following parameter values: 
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The corresponding process window, aerial images and contours on the wafer are given in Fig. 7. Compared 
to the pure "geometrical optimization" in Fig. 6(b) we observe that simultaneous optimization of the source 
further increases the process window. Note that in total we performed roughly 1,000 forward computations for 
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Figure 6. Comparison of two source mask parameter combinations leading to vanishing and large process window, (a), 
(b) Process window; (c),(d) aerial image; (e),(f) developed structure on wafer for intensity levels corresponding to target 
hole width ±10%. Parameters: defocus= — 100 nm, n = r2 — 0.8, b\ = 62 = 45°; (a),(c),(e): w = 192 nm, h = 568 nm, 
d = 104 nm; (b),(d),(f): w = 160 nm, h = 488 nm, d = 40 nm . 
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Figure 7. (a) Optimal process window after source mask optimization; (b) developed structure on wafer at zero defocus 
for intensity levels corresponding to target hole width ±10%. (c),(d) Respective aerial images for zero and — 100 nm 
defocus. Source and mask parameters are given in Equation (9). 

the whole optimization procedure. This is significantly larger than the number of 40 snapshot solutions which 
were computed in the offline phase for construction of the reduced model. 

The presented source mask optimization was performed on a standard desktop PC (computation time of 
the order of one hour, including computation of aerial images and structures on the wafer). Without model 
reduction, the present analysis would have taken over a month. For freeform illumination we expect an even 
larger number of necessary iterations in an optimization procedure such that offline construction of a reduced 
basis model pays off even more. 



5. CONCLUSIONS 



We presented the reduced basis method as a fast and accurate approach to rigorous solutions of electromagnetic 
scattering problems in many-query and real-time contexts. In an offline phase the reduced basis method self- 
adaptively constructs a low dimensional approximation to a given, e.g., geometrically parametrized scattering 



problem. The reduced system is a scattering problem stated on a reduced basis space and can be solved online 
very fast in the order of a second. The foundation for construction of the reduced basis approximation is an 
efficient, rigorous Finite Element solver, e.g., JCMsuitc, 12 which was used in the present work. Since the 
underlying physical model is gencrically incorporated in the reduced model, exponential convergence rates are 
observed, which leads to very low dimensional approximations and relatively small construction times of the 
reduced model. Error estimators furthermore assure the reliability of the reduced basis solutions. 

In a numerical investigation, a reduced model for a geometrically parametrized photomask for cutting hole 
lithography was constructed and incorporated in a source mask optimization project. The constructed reduced 
model could be solved about 10,000 times faster than the full problem. A convergence analysis demonstrated 
that the reduced basis results are in excellent agreement with the exact Finite Element solutions. Furthermore, 
the good performance of the a posteriori error estimator was demonstrated. Estimated and true errors were 
highly correlated. Incorporation of the reduced basis technique into the source mask optimization (SMO) process 
enabled application of combined global and local optimization algorithms for maximization of the process window 
in a seven-dimensional parameter space. The reduced basis offline phase with only 40 forward solutions thereby 
paid-off multiply. After construction of the reduced model the optimization with about 1,000 forward solutions 
could be performed in the order of an hour instead of several weeks, necessary with the full model. 
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